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ABSTRACT 



This thesis presents an analysis of the throughput characteristics for several 
classes of Aloha packet networks. Specifically, the throughput for variable packet 
length Aloha utilizing multiple power levels to induce receiver capture is derived. 
The results are extended to an analysis of a selective-repeat Automatic Repeat 
Request (ARQ) Aloha network. Analytical results are presented which indicate 
a significant increase in throughput for a variable packet network implementing 
a random two power level capture scheme. Further research into the area of the 
near/far effect on Aloha networks is included. Improvements in throughput for 
mobile radio Aloha networks which are subject to the near/far effect are presented. 
Tactical Command, Control, and Communications (C3) s}^stems of the future will 
rely heavily on Aloha ground mobile data networks. The incorporation of power 
capture and the iiear/far effect into future tactical networks will result in improved 
s}'stem analysis, design, and performance. 
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I. INTRODUCTION 



United States Department of Defense (DOD) Tactical Command, Control, and 
Communications (C3) systems of the future will rely heavily on ground mobile 
data networks. These networks must pro\*ide high reliability and connectivity in an 
environment in which enemy counter-C3 assets will attempt to disrupt and deny 
friendly C3 networks. A great deal of interest has been focused on Aloha network 
architectures for future mobile tactical data systems due to their simplicity and 
robustness. 

Research on Aloha has been conducleci since the early 1970s [Kef. l:p. 253]. Most 
of this work. liowe\'er, has utilized idealized channel models. While this has pro- 
duced a wealth of information about the general characteristics of Aloha networks, 
it remains critical that real world conditions be taken into account. Factors which 
afi’ect Aloha channels in the real tactical environment include fading, Near/Far ef- 
fect. pow('r ca]^tur<* (both natural and induced), and jamming. This thesis presents 
analysis of Near/Far efb'ct and ])ower capture' on channels which employ the Aloha 
network discipline'. 

Chapters 11 and III concentrate on the effects of induced power capture when 
employed as part of a \'ariable ]>acket Aloha network. Chapter I\' presents an anal- 
ysis of the ramifications of Near/Far effect on three different Aloha configurations. 
This work is geared specifically toward ground mobile data network applications. 
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II. POWER CAPTURE ALOHA 



Aloha packet networks have been utilized for some time in data communications 
systems which are composed of many uncoordinated users, each with low data rate 
requirements. Recent emphasis on military and commercial ground mobile data 
networks will increase the number of systems which employ variations of the Aloha 
protocol as a channel discipline. Although many Aloha systems have been effec- 
tively employed to date, several variations on the basic Aloha scheme still remain 
unstudied. The preponderance of work on Aloha network analysis has focused on an 
idealized theoretical radio channel, thereby neglecting the effect of real world chan- 
nel conditions on system performance. It is particularly important to incorporate 
channel characteristics into Aloha analysis because, in contradiction to intuition, 
the overall throughput of an Alolia channel can actual!}' increase under certain con- 
trived or natural channel conditions. A specihc example of this is the case of power 
capture Aloha. 

Another area of neglect in terms of Aloha network analysis is that of variable 
packet Alolia. Although model results have shown that fixed packet and slotted 
Aloha both exhibit better throughput characteristics than variable packet, the real- 
ity of bit stuffing logic in fact reduces the actual throughput of these disciplines. In 
contrast, variable packet Aloha can be shown to be a simple and efficient network 
discipline. 

This chapter introduces and develops the throughput characteristics of a variable 
packet Aloha network which incorporates the power capture phenomenon. For this 
analysis, only created capture is discussed, although the results can easily be adapted 
to accommodate natural capture due to fading. 
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Created capture effects occur when groups of users are assigned different transmit- 
ter power levels in an effort to create priority classes, or when all users dynamically 
select a random transmit power in order to increase channel throughput. Previ- 
ous analysis of power capture in slotted and fixed packet networks has shown a 
dramatic increase in channel throughput over the idealized channel models [Ref. 2]. 
This chapter investigates the effect of created power capture with two random power 
levels on an Aloha network using variable length packets. 

A. SYSTEM DYNAMICS 

In this analysis, it is assumed that there is an infinite user population and that 
the channel traffic rate is Poisson with parameter g packets per second. Assume 
newly created packets are of variable length x (bits) and are independent!}’ and 
identically distributed according to the probability density function n(a’) with mean 
X. Whthin this model, all channel packets are indistinguishable. That is to say that 
newly arrived packets and retransmitted packets have their length x redrawn afresh 
from a (a*) for every packet . This is an analytical device known as the independence 
assumption. 

The probability that n packets interfere (overlap) with an arbitrary tagged packet 
during its transmission is given by 

/>,(„) = + ,2 1 ) 

where x is the tagged packet length, u represents the random length of a preceding 
packet, and n is the number of interfering packets [Ref. 3]. For the case of n — 0, 
this can be interpreted as the probability that the tagged packet encounters no 
overlap from pi'eceding packets and the next interarrival is larger than x. Given the 
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case where no interferers can be tolerated by a tagged packet, equation (2.1) can be 
viewed as the probability of successful tagged packet transmission. 

This analysis, however, considers the case when a receiver can be captured by a 
tagged packet in the presence of interfering packets if the tagged packet signal power 
exceeds the joint signal power of all interfering packets by some specific capture 
threshold, 70. If at any time during the tagged packet interval, the ratio of tagged 
packet power to joint interfering power falls below 70, the tagged packet is considered 
to be destroyed. 

Unlike slotted Aloha where the number of interferers is constant over the entire 
tagged packet transmission interval (Figure 2.1(a)), or fixed packet Aloha where the 
number of interferei's is a stochastic process made up of early and late interfering 
packets (Figure 2.1(b)). the number of interferers present in a variable packet Aloha 
channel represents a stochastic process which is driven by four classes of interferers. 
Figure 2.1(c) depicts the four possible interferers types. The terms early and late 
refer to the begin and end time relationship between an interfere!' and the tagged 
packet. For instance, the interfere]' identified as Early-Early in Figure 2.1(c) cor- 
responds to an\' iiilerferer which begins before the tagged packet (early) and ends 
before the tagged packet (early). 

The dynamics of the interfering stochastic process is as follows. We note that the 
arrival and departure of n interfering packets partitions the tagged packet interval 
into several non- overlapping intervals of random length. The exact number of these 
intervals, gi\’en a known number of interferers and their corresponding interfere!' 
type, is found by the following relationship. 

Total Intervals = (number of Late-Late packets) 

+ (number of Early-Early packets) 
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Figure 2.1: Aloha Interfering Packets 
+ 2 X (number of Late-Early packets) + 1 

Relating this equation to Figure 2.1(c), it is evident, for instance, that a single Early- 
Early or Late-Late packet would partition the tagged packet into l-fl intervals. One 
interval occurs during the absence of the interferer; the other during the presence 
of an interferer. In a different case, a single Late- Early would partition the tagged 
packet into 2 -f 1 intervals. Two of the partitions would have no interferers present; 
the other partition would occur during the period of the interfering packet. The 
Early-Late packet covers the entire duration of the tagged packet and therefore 
does not provide any partitioning. 

Each interval has associated with it a total number of interferers /(/,), where t, is 
the start time of the Ftli partition during the tagged packet interval. Each uniciue 
sequence of total interferers per partition (/(/,) for all i) shall be called a realization. 
Figure 2.2 depicts a specific realization of interfering packets for a case where v = .3. 




Figure 2.2: An Example of an Interference Realization 

For this specific realization, the tagged packet total intervals is 5. The total number 
(/(/,)) of interfeiers present during each interval is indicated at the bottom of the 
figure. 

B. MAXIMUM NUMBER OF INTERFERERS 

A computer program ( A])j)cndix A) was written to generate a list of every i)ossibIe 
interfering realization for a given ??. From this list of all possible realizations, the 
number of realizations {Cj) in which the maximum number of interfeiers present 
at an}^ time during the tagged packet interval does not exceed some number j was 
tallied. Table 2.1 gives the results for n = 1,2, 3,4,5. 

An example of how one would generate a Cj{n) value is as follows. Consider 
the simple case of n — 1 . In this case, only four realizations are possible, namely 
the interfering jiacket can be an Early- Early, Early- Late, Late- Early, or Late- Late. 
Continuing with this simi)le examj^le, it can be seen that given there is 1 known 
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n 


Ci(n) 


C2{n) 


C3{n) 


C-i(n) 


Cs{n) 


It 




1 


4 


0 


0 


0 


0 


4 


0.25 


2 


4 


10 


0 


0 


0 


14 


0.0714 


3 


4 


24 


20 


0 


0 


48 


0.0208 


4 


4 


50 


75 


35 


0 


1C4 


O.OOGl 


5 


4 


100 


225 


17C 


5G 


5C1 


0.0018 



TABLE 2.1: Maximum Number of Interfering Packets Given n 

interfercr, there is no realization where the niaxiinuin number of siiiiiiltaiieous in- 
terferers during the entire duration of the realization is less than or equal to j = 0. 
In fact, in this example, the maximum number of simultaneous interferers for each 
of the four possible realizations is 1. Therefore, the value of Ci{n) for n 1 is 4 
as seen in Table 2.1. For (he case of n = 2, there are 14 combinations of 4 packet 
types taken 2 at a time. These 14 pairs of interfering packets generate 14 unique 
realizations. Of these 14 realizations, 4 have a maximum simultaneous number of 
interferers of 1, and 10 have a maximum of 2 interferers. The program in Ap])endix 
A is used to automate the process of generating combinations of packets, construct- 
ing the resulting realizations, and determining the maximum number of interferers 
per realization. Table 2.1 is the result of the program. 

Each Cj{n) gives the number of realizations in which the inaxiinum number of 
interferers at any time equals j, given the total number of interfering packets is 
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known to be n. The function Ij is simply the total number of all possible realizations 
given n packets are known to interfere. The function /(n) is the inverse of Ij and 
is simply used as a shorthand notation for 1//t. 

C. VARIABLE PACKET ALOHA WITH RANDOM POWER LEV- 
ELS 

Data packets arrive at the receiver after being spatially attenuated. Each packet 
is transmitted with one of two normalized power levels given by the set D = {1, M}. 
Here, M is some integer multiple of the threshold 70 times the lower power level 
P = 1 (normalized), that is M — V70 where is an integer greater than 1. Each 
user selects a transmit power level for each individual packet from the set D with 
equal probability | in order to avoid any class prioritization among user groups. 
The higher power level M is chosen according to the relation 

(AT l)7o > M > V70 (2.2) 

where A > 1 is an integer and 70 is the power capture threshold of the receiver. 

The tagged packet, arriving at the receiver with power Pt 6 {l,M} ma}’ capture 
the receiver given a realization of n interfering packets if and only if 

HU) 

Pt > 7omax[^ Pj] (2.3) 

j=i 

for all intervals (i’s) in the realization where Pj is the power level of packet j. 

For a tagged packet to capture the receiver, two events must occur. First, the 
tagged packet must be of power M while none of the interferers are of power M. 
Call this event A, then 



Pr{A\n} 



1 

2("+i) 



(2.4) 
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The second condition for capture is that the tagged packet must satisfy (2.3) where 
Pi = M and all Pj = \ . Given a specific realization, this equates to 

N > max{l{ti)) (2.5) 

The value N is therefore the maximum number of interferers that can be tolerated 
at any instant of time during the tagged packets transmission interval. If event B 
denotes the tagged packet power satisfying (2.5), then 



Pr{B\rt] = f{n) ^ C,{n) (2.6) 

j=0 

where Cj(?i) is the number of realizations, given interfering packets, in which the 
maximum number of simultaneous interferers equals j. Table 2.1 provides a list of 
c,{ n)'s and the inverse /(n) of the total number of possible realizations given n 
interfering packets. 

By using equations (2.4 ) and (2.6) the conditional capture probability Pr{capture\n } 
is obtained. 



Pr{capiurt\n,n = 0] = 1 (2.7) 

J^r{copiurc\n.n>l} = Pr{A\n) Pi'{B\n) (2.S) 

= (2.9) 

J=0 

The overall probabilit\’ of receiver capture by a tagged packet is found by mul- 
tiplying the conditional probability in equation (2.9) by the probability that there 
are 7? packets in the channel, and summing over all n. Equation (2.1) provides the 
probability of n packets, therefore, the probability of receiver capture b\' a tagged 
packet which shall be called the probability of success is 
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Pr{succtss] = (2.10) 

n=l ”• j=0 

The conditional channel departure rate in the interval a: + n is given by Reference 3 
to be 



(f>{x,u) = g{x + u)Pr {success} (2-11) 

The channel throughput S, in packets per mean packet time, is related to the channel 
departure rate by 



roc- roo 2 

S= / -(p{x^u)a{x)a{n)(lxdu (2.12) 

Jo Jo 2 

Define G to be the attempted packet rate in packets per mean packet time (i.e. 
G = gx). Appendix B contains the detailed solution of the above equation. The 
resulting channel throughput is found to be the following: 



N 



s = 



+ t ^^|^/('0G-'(G + !)-'”«> gc,(„) (2.13) 



n=l ^ j=0 

Figure 2.3 shows the throughput of variable packet Aloha with random two power 
signal for values of A" — 0,1, 2, 3. Packet lengths are exponentially distributed and, for 
computational purposes, n is taken to be between 1 and 5. The curve corresponding 
to N = 0 is for conventional variable packet Aloha. 
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VARIABLE PACKET ALOHA UITH POWER CAPTURE 




Figure 2.3: Throughput of Variable Packet Aloha With Power Capture 



11 



III. SELECTIVE-REPEAT ALOHA WITH 

CAPTURE 



To continue the study of variable packet Aloha, the results of the previous chapter 
are now applied to a more advanced Aloha network which incorporates selective- 
repeat as a data link control. 

The specific architecture discussed here is a network where remote terminals com- 
municate with a hub station via an Aloha channel in the query-response mode. The 
feedback hub-to-remote link is an asychronous time-division-multiplexed (ATDM) 
channel. Each remote terminal receives the same data stream and selects applicable 
packets by searching for its own address. 

Assume that packet length of newly arrived and retransmitted packets in the 
Aloha channel are independent and identically distributed with probability density 
function a(x) and with mean x. In this model, packets are indistinguishable. Both 
newly arrived packets and retransmissions can be pictured as having their length 
redrawn afresh from a{x) prior to each transmission. 

Packets in the ATDM channel are independent and identically distributed with 
density function b{y) and mean y. 

Variable packets on the Aloha channel are broken up into ^ numbered and con- 
stant length minipackets prior to transmission. When packet damage occurs, only 
the specific minipackets which are unreadable are retransmitted. Let D{z)^z =1,2,.. 
be the distribution of the minipackets. On the average, each large packet is seg- 
mented into ^ minipackets. 
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Let the channel attempted packet rate be g packets per second and 'zg the channel 
attempted rate in minipackets per second. With m being the fixed length of a 
minipacket, the attempted channel rate in minipackets per minipacket time is mzg. 

The conditional probability that a minipacket encountered no collisions is the 
probability of no overlap from preceding packets and that the next interarrival is 
larger than the minipacket size m. Thus 

P,{z) = (3.1) 

In a capture environment, however, Ps must take into account the times when the 
tagged packet captures the receiver in spite of the presence of n interferers. There- 
fore, using equation (2.1) from the previous chapter, the probability of minipacket 
success is 



P.(z) = + T y c^(„)] (3,2) 

n! 



7t = l j — 0 

The conditional probabilit}’ of successful minipacket transmission is 



q{z,y) = P 5 ( 2:)[1 -Pi(m)][l -p 2 (y)] (3.3) 

where Pi(m) is the minipacket error probability due to errors in the Aloha channel 
and P 2 {y) is the error rate in the ATDM channel. The conditional minipacket 
departure rate in the interval of r + 1 minipackets is 



4>m(:-y) = mg{z + \)<l{z,y) 



(3.4) 



The average minipacket throughput is giv^en by the following equation: 



z 1 Jo 



4>m{~zy)dy 



(3.5) 
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Assume that each packet consists of a geometrically distributed number of fixed 
length minipackets, then the distribution D{z)^ which is the probability that a packet 
consists of 2 : minipackets, is defined as follows. 



C(2)=(l-,<rV (3.6) 

where // is the probability that a minipacket is generated. Because 2 = D{z) 

can be written as 



= hi - h”" 



(3,7) 



Appendix C provides the detailed solution for the throughput equation. The 
resulting rather length}’ set of equations is as follows: 



G 



= A + - — <7(1 -pi{m)) h{n)w ‘^[F{n)R{n) - /?] (3.8) 

“r ->■ 1 



n=l 



A = 



aGv^{a — r’)(^ ~ 
z(2 T l)u;^(l — vY 



(3.9) 



F{n) = (1 -u)-('+’)(t + 1)! 



(3.10) 



t r— 1 






r = l /rz=0 



(3.11) 



E{k,r) = 



jr - ky 

k\{t-k + l)\ 



(3.12) 



G ■ 7nzg 



(3.13) 
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(3.14) 



h{n) = 



G" 



Ec,(") 



j=0 



= (1 - 


(3.15) 


«; = (!-=) 

z 


(3.16) 


k = z 1 


(3.17) 


i = n 1 


(3.18) 



Figure 3.1 shows the throughput for the \^ariable packet Aloha with selective^ 
repeat data link control. The parameters are: P\{rn) = 0, = 0, x and y are 

both exponentially distributed, Y = y, and z = S minipackets. 
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SELECTIVE REPEAT UITH POUER CAPTURE 




Figure 3.1: Throughput of Selective- Repeat Aloha with Power Capture 
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IV. NEAR/FAR EFFECT ON ALOHA 

MODELS 



This final analysis chapter deals with a network phenomenon which is specific to 
mobile packet radio networks that do not incorporate adaptive control of terminal 
transmit power. The near/far effect is a condition which occurs when packets arrive 
at a receiver with different mean power due to differences in transmission distances. 
Under these conditions, a tagged packet may be able to capture the receiver, in spite 
of interfering packets, if the tagged packet power is sufficiently greater than the joint 
power of the interfering packets. Because each packet is spatially attenuated due 
to the iiear/far effect, all packets have a finite probability of capturing the receiver. 
The throughput characteristics of various Aloha network configurations which are 
subject to the near/far effect are derived and presented. 

A- CELL MODEL 

Idle mean power tr of a packet at a distance from the transmitter is of the 
general form 



w = cr ^ (4.1 ) 

In the event of ground-wave propagation without shadowing. 

c = Pt^Gt^rHIhI (4.2) 

where Fj,, Gj,. and Hji are the transmit power, antenna gain, and antenna height 
(above ground), respectively, of the transmitting mobile terminal. The values Gr 
and Hfi are the gain and heiglit above ground of the base station antenna. The 
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exponent a gives the attenuation law for the channel considered (2 < a < 5). In 



the event of UHF propagation, a typical value is a = 4 [Ref. 3:p. 264]. 

Assuming identical mobile terminals and omnidirectional antennas c is normal- 
ized to unity for all users without loss of generality since receiver capture power is 
determined by the ratio of signal powers [Ref. 3:p. 264]. Additionally, the value of 
Q in this analysis is taken as 4. 

In this model, signal power depends on the spatial distribution of users in the 
network. Users create a trafhc density g{r) within a circular cell around a base 
station receiving node. The density ^( 7 ^) can be thought of as the normalized packet 
traflic per unit area at a distance 7’ from the base receiver. The packet generation 
rate in a given area depends on the distance to the base station r and is independent 
of direction. Given this cell model, the total traffic rate offered in the network can 
be described as 




(4.3) 



Figure 4.1 depicts the general cell model used in this analysis. 



The spatial distribution of packet generation is found to be the probability a 
packet is generated within a distance R. That is to say 




(4.4) 



Differentiating, the probability density function for packet generation is found to 



be 




( 4 . 5 ) 



IS 



ALOHA NETWORK 




Figure 4.1: Cell Model for Near/Far Effect Analysis 

As a means of siinplificatioii, users are assumed to be uniformly distributed in a 
circle of radius 'vitli the base iecei\'iug station at tlu' center of the cell as shown 
in Figure 4.2. As seen in Figure 4.2, 5 ( 7 ') is normalized to for ()<?'< 

From this, the equations for the total network traffic and packet generation reduce 
to 



. , , 27r7‘ 

fn{r) = (4.7) 

for 0 < r < 

A standard change of variable operation is performed on equations (4.1) and (4.6) 
to define a probability density function for the packet power w. 



19 



g(r) 



G 

n r 



2 



max 



0 



-I r 
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Figure 4.2: Distribution of Users Within the Cell 

fw{w) = 

for < u; < oo. 

At this point, the network cell model has been described with the distribution 
of packet generation given by equation (4.6) and the packet power density function 
given by equation (4.7). The near/far effect is now incorporated into the model. 

B. NEAR/FAR EFFECT 

In a phenomenon which is similar in many respects to the power capture case 
described earlier, the near/far capture occurs when the ratio of a tagged packets 
power to the joint interfering power exceeds some receiver capture threshold. Again, 
the number of interferers during the tagged packets transmission is a stochastic 
process. The capture ratio is described as 
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(4.9) 



Z = 



X 

Y 



where X is the packet power of the tagged packet and Y is the joint interfering 
power of some number n interfering packets. As in previous discussions, the receiver 
is captured if the ratio described by Z exceeds the threshold given by the parameter 



7o- 

The distribution of A' is simply the single packet power density function given 
in (4.7). The joint interfering power density is dependent on the number of 
interferers and is given as 



fy(y\n) = [/ir(u’)]"® (4T0) 

where iiQ denotes the //-fold convolution. For the cases of 7/ = 0 and 77 = 1, no 
convolutions are required to determine fy{y). For higher numbers of interferers, a 
computer based convolution package is used to generate specific values for fy{y\7i)- 
Appendix D provides details on the anah^tical result for the case of 77 = 1 and the 
numerical results for u = 2.3. 

Utilizing the density functions for and Y, the distribution of Z can be found 
using the following formula. 



TOO 

fzi = \n)= yfx{y~)fY{y\n)dy (4.11) 

Jo 

The capture probability, that is the probability that Z exceeds the receiver thresh- 
old 7 o given n interferers. is found by integrating fz{^\n) from 70 to 00 . 

Table 4.1 presents the probability of capture given n interferers for the range 
0 < n < 4. Although the case of /; = 1 is in a closed form with parameter 70, all 
subsequent capture probabilities assume a value of 70 = 3. 
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n 


Pr{caplurc\n } 


0 


] 


J 


•pu 


2 


0.057 


3 


0.022 


4 


O.OJO 



TABLE 4.1: Capture Probabilities 

The limited list of capture probabilities listed in Table 4.1 is sufficient to provide 
excellent resolution of the throughput characteristics for the Aloha models discussed 
in the following sections. The probability of capture for cases where there are greater 
than 4 interfering packets is so small as to be considered insignificant. 

C. ALOHA NETWORK MODELS 

ddiree common Aloha models are analyzed with respect to the improvement in 
overall network throughput gained from the near/far effect. In each case, the general 
definition of throughput is 



S = GPr {success] 

which in a capture environment can be viewed as 



(4.12) 
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(4.13) 



5 = G ^ Pr{capture\n] Pr{n} 

n=0 

The following sections use this approach to determine the network throughput for 
slotted and fixed packet Aloha networks. Variable packet Aloha is analyzed in the 
same general manner as was done in Chapter I for power capture. 

1. Slotted Aloha 

For slotted Aloha, the probability that n packets are generated during a 
given packet time is given by the Poisson Distribution [Ref. l:p. 255]: 

Pr{n] = (4.14) 

Therefore, the general form for the throughput in a slotted Aloha network with the 
near/far effect is 



<x G^e~^ 

S = F7'{capture\n] ; — (4-15) 

74=0 

Using the Pr{captiire\n] values in lable 4.1. the throughput S is calculated and is 
shown in Figure 4.3. 

2. Fixed Packet Aloha 

The number of packets generated during a given packet time in a fixed 
packet network is given by [Ref. 4:p. 11]: 

Pr{7r} = (4 jg) 

n\ 

resulting in a throughput of 

00 

.S’ = G y] pr {capture\n 

71 = 0 






n\ 



( 4 . 17 ) 
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SLOTTED ALOHA WITH NEAR/FAR EFFECT 




Figure 4.3: Slotted Aloha With Near/Far Efiect 
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FIXED PACKET ALOHA UITH NEAR/FAR EFFECT 




Figure 4.4: Fixed Packet Aloiia With Near/Far EfTect 

Again, using tlie miinorical results for the probability of capture given in Table 
4.], the tl)rougb|)ut is cak iilated and is de]>icted in Figure 4.4. 

3. Variable Packet Aloha 

This analysis is very similar to the work j)rescntcd in Chapter II and 
therelore the variable names used here are the same. First, the jrrobability of success 
for a tagged packet is defined as 



P,{x,ii) = Pr{capture\n)e (4.18) 

n=o ”• 



2.5 



For this analysis, perfect channel conditions are assumed in order to simplify the 
calculations. The conditional probability of a successful transmission in this case is 
seen to be 



q{x,u) = Ps{x,u) (4.19) 

The conditional channel departure rate in the interval a; + u is 

4>{x, u) = g{x + u)q{x + u) (4.20) 

From this, the channel throughput can be calculated from the following equation: 

S = - / (j){x^u)a{x)a{u)dxdu (4.21) 

2 Jo Jo 

Take the case where a{x) and a{u) are exponentially distributed, the above equa- 
tion simplifies to 

^ {G+lf ^ T ” ^ ~ T r { capture | n } ^ ( G + l)"("+3) (4.22) 

Using the numerical results for Pr{capture\n} given in Table 2.1, the throughput 
for variable j^acket Aloha with Near/Far effect is shown in Figure 4.5. 
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VARIABLE PACKET ALOHA WITH NEAR/FAR EFFECT 




Figure 4.5: Throughput for Variable Packet Aloha With Near/Far Effect 
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V. CONCLUSIONS 



This thesis has demonstrated the effects of induced two level power capture on 
variable packet Aloha networks and the effects of the near/far phenomenon on var- 
ious Aloha configurations. In all cases, significant improvements in overall network 
throughput were achieved. 

The development of highly reliable Tactical data systems of the future depend 
on Aloha as a multiaccess scheme for ground mobile terminals. It is critical that 
these networks be allowed to achieve their maximum potential. Given the increased 
reliance on tactical data conununications, Aloha networks which provide the maxi- 
mum user throughput possible, even in a jamming environment, must be developed. 

The increases in overall network throughput presented in this thesis indicated 
the improvements in performance which can be realized when power capture and 
the near/far effect are incorporated. As seen in Chapter 2, induced power capture 
increased the network throughput of a variable packet Aloha network from 0.14S 
packets per mean packet time in the idealized case to 0.184 in a capture environment. 
Utilizing the selective- repeat data link control with power capture variable packet 
Aloha increased throughput from a maximum of 0.256 to about 0.354. In the final 
analysis chapter, improvements in throughput for slotted, fixed packet and variable 
packet Aloha which incorporated the near/far effect were presented. 

Continued work in the area of Aloha network analysis will be needed in order 
to develop the high capacity, high reliability tactical data networks required in the 
future. Incorporating power capture and the near/far effect into future designs will 
result in better systems anah’sis and design for these tactical mobile networks. 



APPENDIX A 



{ } 

{ PROGRAM CJ } 

{ } 

{ Joe McCartin, January 1989 } 

{ } 

{ Turbo Pascal - IBM PC } 

{ } 

{ Description; } 

{ } 

{ This prograjn is designed to generate all } 

{ realizations of interferers given n packets are } 

{ known to interfere. Using the realizations, a } 

{ table of Cj (n) values is compiled. This progreim } 

{ is geared specifically toward the purposes of this } 

{ thesis and is not a general purpose program. } 

{ > 

program cj (filevar, output); 

const 

max_intervals = 11; 

matrix_length = 1024; 

cj .matrix_dim = 10; 

type 

realization_vector = array [l . .max_intervals] of 

integer ; 

realizat ion_matrix = array [l . .matrix.length] of 

realization^vector ; 

cj^vector = array [l . . cj _matrix_dim] of integer; 

cj.array = array [l . . cj_matrix_dim] of cj_vector; 

index = array [1..4] of integer; 

comb_matrix = array [1..220] of index; 

var 

realizat ion_x , 

realizat lon^y : realizat ion^matrix ; 

vector : comb_matr ix ; 

vector_length : integer; 

act ive^matrix , 

act ive_mat rix.length : integer ; 



cj_matrix : cj^array ; 



init^interf erers , 
max.interf erers , 
min_interf erers , 
num.intervals , 
final.interf erers : integer ; 

n, 

i, 

j. 

t , 
h, 

runs : integer; 

file_name:string[lO] ; 
f ilevar : text ; 
outf ilevar : text ; 

ans : char ; 

archive : boolean ; 



{ } 

{ INITIALIZE.CJ.MATRIX } 

{ } 

{ Tabulated values of Cj^s are kept in a square matrix } 
{ defined by the dimension cj _matrix_dim . This proc } 
{ initialize this matrix which is called CJ^MATRIX. } 

{ } 



{^:4c:*c**:f::f:***********Jf:*******5f:3*c***********itc******:*c*:f:******} 

procedure INITIALIZE_C J^MATRIX ; 
var 

1 , 

j : integer; 
begin 

for i:= 1 to cj_matrix_dim do 
begin 

for j:= 1 to cj _matrix_dim do 
begin 

cj.matrix [i,j] := 0; 

end; 

end ; 

end ; 



{ } 

{ INIT.REALIZATION MATRIX } 

{ } 

i In this program, realizations are generated one } 
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{ interval at a time until the entire realization is } 
{ constructed. For each interval, a matrix of all } 
{ known realizations is read. Additions to these } 
{ partial realizations are then placed in a different } 
{ matrix. This proc initialize both of these } 

{ matrices. } 

{ } 



procedure INIT_REALIZATION_MATRIX ; 
var 

i, 

j : integer; 
begin 

writeln (INITIALIZING REALIZATION MATRICES '); 
writeln; 

for i := 1 to matrix.length do 
begin 

for j := 1 to max_intervals do 
begin 

realization.x [i , j] := 0; 
realization_y [i , j] := 0; 
end; 

end ; 

writeln (INITIALIZATION COMPLETE 0; 
writeln; 
end ; 



{ 

{ C0PY_X_T0_Y 

{ 

{ This procedure copies the contents of realization 
{ matrix X into realization matrix Y. 

{ 

procedure C0PY.X_T0_Y ( x ,y : integer) ; 
var 

i : integer; 
begin 

for i := 1 to max_intervals do 

realizat ion^y [y , i] := realization_x[x, i] ; 

end ; 
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{ } 

{ COPY Y TO X } 

{ > 

{ This procedure copies the contents of realization } 

{ matrix Y into realization matrix X } 

{ > 



procedure C0PY_Y.T0_X (y,x : integer) ; 
var 

i : integer ; 
begin 

for i := 1 to max. intervals do 

realization.x [x, i] := realizat ion.y [y , i] ; 

end; 






{ } 

{ PRINT MATRIX } 

{ } 

{ Print.matrix prints out the most recently updated } 

{ realization matrix (X or Y) . If the flag MATRIX is } 

{ set to 1, then the REALIZATION. X was the last } 

{ updated and is subsequently printed. If the flag } 

{ is not 1, the Y matrix is printed. If the } 

{ diagnostic flag ARCHIVE is true, output goes both } 

{ to the screen and to the data file OUTFILEVAR. } 

{ Else, output is to screen only. } 

{ } 



-^:f:;if::4c:4c3f:5tc:4c:+:5*c:4:;tc5f::4c:4c3*cj*e:4c5f::*c:4:5f:5f::f:5f:3f::+:3f:3tc;f:jf:;f::4c;tc:4c:4c:4c*;f:3f:*:4cjf:*;f:5f:* 

procedure PRINT.MATRIX (matrix, length, width : integer) ; 
var 

i, 

j : integer; 
begin 

for i := 1 to length do 
begin 

writeln; 
if archive then 

writeln (outf ilevar) ; 
for j := 1 to width do 
begin 

if matrix = 1 then 
begin 

write (realization.x [i , j] ) ; 
if archive then 
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write (outfilevar, 

realization_x[i, j] ) ; 

end 

else 

begin 

write (realization_y [i , j] ) ; 
if archive then 

write (outfilevar, 

realization.y [i , j] ) ; 

end; 

end; 

end ; 

writeln; 
if archive then 

writeln (outfilevar); 

end; 



■(********♦**+**+****♦*****♦*****♦♦**♦**♦****♦♦♦♦*♦**♦**} 



{ } 

{ GENERATE.REALIZATIONS } 

{ } 

{ This procedure produces a matrix containing all } 

{ possible realizations given an initial number of } 



{ interferers, a maximum number of interferers, and } 
{ the total number of intervals. A tree method is } 
{ employed to build up a realization one interval at } 
{ a time. The first interval contains the initial } 
{ number of interferers. This partial realization } 



{ is placed in a realization matrix. On the next } 
{ pass, the partial realization is read and the next } 
{ interval is added to it. These realizations are } 
{ then placed into another realization matrix. This } 
{ process continues until all intervals have been } 
{ filled in for all realizations. } 
{ } 



procedure GENERATE^REALIZATIONS (init.interf erers , 

max. interferers , 

num.intervals , 
min. interferers : integer; 

var active.matrix : integer 
var active. matrix. length 
: integer) ; 



temp, 

i, 

J. 

k, 
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var 



x_list_length, 
y.list^length, 
interval^count : integer; 

begin 

writeln (^GENERATING REALIZATIONS FOR \m 2 LX_interf erers , ^ 
MAX INTERFERERS AND \num.intervals / INTERVALS^; 

{ writeln (init^interf erers, ' \max_interf erers , ^ 

^ ,nuin. intervals, ^ ^ ,roin_interf erers) ;} 
writeln; 

if archive then 
begin 

writeln (outf ilevar) ; 
writeln (outf ilevar , 

^REALIZATIONS ************** ^ ) ; 

end ; 

realization_x [l , l] : = init_int erf erers ; 
x_list_length := 1; 
interval.count := 1; 
active_matrix := 1; 
active_matrix_length := 1; 

while interval_count < num.intervals do 
begin 

y_list_length := 0; 

{ writeln (^within x to y loopO;} 

for i := 1 to x_list_length do 
begin 

if realization_x [i , interval_count] < 
max_interf erers then 
begin 

y_list_length := y_list_length + 1; 
copy_x_to_y (i ,y_list_length) ; 

realization_y [y_list_length , interval_count+l] 

: = realization^x [i , interval_count] +1 ; 
end; 

if realization.x [i , interval_count] > 
min.interf erers then 
begin 

y_list_length := y_list_length + 1; 
copy_x_to_y(i,y_list_length) ; 

realizat ion.y [y^list^length , interval_count+l] 

: = realizat ion_x [i , interval .count] -1 ; 
end; 



end; 

interval.count := interval.count + 1; 
active.matrix := 2; 

active_matrix_length : = y_list_length ; 
x_list_length := 0; 

if interval. count < num.intervals then 
begin 

for k := 1 to y. list. length do 
begin 

{ writeln within y to x loopO;} 

if realizat ion.y [k , interval. count] 

< max.interf erers then 

begin 

X. list. length := x.list.length 
+ 1; 

copy. y_to.x(k, x.list.length) ; 

realizat ion. X [x.list.length , interval. count +l] 

:= realization.y[k, interval. count] +1 ; 
end ; 

if realization. y [k, interval.count] 

> min. interf erers then 

begin 

x.list.length := x.list.length 

+ i ; 

copy.y.to.x(k, x.list.length) ; 

realizat ion. X [x.list.length , interval. count+l] 

:= realization. y [k, interval.count]-! ; 
end ; 

end; 

interval.count :=> interval.count + 1; 
active. matrix := 1; 

active. matrix. length := x.list.length ; 
end; 



end ; 

writeln (^There Were \active. matrix. length , ^ 
Realizations ' ) ; 
if archive then 
begin 

writeln (outfilevar /There Were 
^ , active. matrix. length, ^ Realizations/ ; end; 

if active. matrix = 1 then 

print. matrix (l , act ive. matrix. length, 
interval.count) 

else 

print. matrix (2 , act ive. matrix. length, 
interval.count) ; 
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end ; 



•^:tc3|c:f4(:fc:f:4c:4c;itc:tc:tc:f4(]1c]1c:4c]tc:4(:tc]4e:4(3tc3lc:4c:tc]1c:4(3tc]1c:4c]tc:4c:fc:tc:f3tc3tc;1c:fc:fc:f:f(3fc:f:tc:tc]tc:4c:1c3lc:4c]1c:4c:4c^ 

{ > 

{ FILTER } 

{ } 

{ This procedure accepts a realization matrix which } 
{ contains all possible realizations given a specific } 
{ combination of interferers. Each realization in } 
{ the matrix is checked for the invalid condition of } 
{ overcount. An overcount occurs when a realization } 
{ has indicated that a certain number of interferers } 
{ has departed but yet somewhere in the realization } 
{ the interferer count goes beyond the count of the } 
{ remaining interferers. Realizations which are } 

{ filtered out are discarded. All good realizations } 
{ are placed into the currently inactive realization } 
{ matrix. } 

{ } 



procedure FILTER (var length: integer ; intervals, 
interferers : integer; 

var active_matrix : integer) ; 



var 

i, 

1 > 

P. 

max , 

num.f inal.interf erers : integer; 

invalid_realizat ion : boolean ; 

begin 

writeln; 

writeln ('FILTERING THE REALIZATION MATRIX '); 
writeln; 



j := 1; 

if active.matrix = 1 then 
begin 

active. matrix := 2; 
for i := 1 to length do 
begin 

invalid. realization := false; 
num.final. interferers : = 

realization.x [i , intervals] ; 
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if num_f inal.interf erers = interferers then 
begin 

if intervals > 2 then 
begin 

mcLX := max_interf erers ; 
for p := 1 to (intervals - 1) do 
begin 

if realization.x [i ,p] > 

realization^x [i ,p+l] then 
max := majc - 1; 
if realization_x [i ,p+l] > max 
then 

invalid.realizat ion 
:= true; 

end ; 

end ; 

if invalid.realization = false then 
begin 

copy_x_to_y 

{ print_matrix(2, j) ;} 

j := j + 1; 

end ; 

end; 

end; 

end 

else 

begin 

act ive_matrix := 1; 
for i := 1 to length do 
begin 

invalid^realization := false; 
num_final_interf erers : = 

realization_y[i , intervals] ; 
if num_final_interf erers = interferers then 
begin 

if intervals > 2 then 
begin 

max := max_interf erers ; 
for p := 1 to (intervals - 1) do 
begin 

if realization_y [i ,p] > 
realization_y [i ,p+l] then 
max : = majc - 1 ; 
if realization_y [i ,p+l] > max 
then 

invalid, realization 
;= true; 

end ; 

end ; 

if invalid_realization = false then 
begin 
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{ 



copy_y_to_x 
print_matrix 
j := j + 1; 
end; 

end; 

end; 

end; 

writeln; 

writeln (^THERE WERE VALID REALIZATIONS FOUNDS; 

writeln; 
if archive then 
begin 

writeln (outf ilevar) ; 

writeln (outf ilevar / FILTER ******^******:^**:^*** ^ ) 
writeln (outf ilevar /There Were /j-1/ Valid 
Realizations ^ ) ; 
end; 

length := j”l; 

print .matrix (active.matrix, length, intervals) ; 

writeln; 

if archive then 

writeln (outf ilevar) ; 

end ; 



{ } 

{ GET.Cj } 

{ > 

{ GET_Cj reads in all valid realizations individually } 

{ and determines the maximum number of interferers } 

{ in each realization. A tally is kept of how many } 

{ realizations have a given maximum number of } 

{ interferers. } 

{ } 

-[J*c*3f:3tc:4c3tc3tc5^:*c:f:5f:3*c:tc5*c^:4c3tc:^*:tc:+::!f:3tc3tc)tcj*c3*c:4c:*c:+::4c:4c3*c;4c3tc3tc3tc3*c)tc:4c3*c3tc3t::tc:tc:4c:tc;tc3*c3tc:+::tc;f:3f:)- 



procedure get.cj (active. matrix , act ive.matrix. length , 

n: integer); 



var 

i, 
j > 

max : integer ; 
begin 

writeln (^FINDING Cjs FOR A REALIZATION WITH /n,^ 
PACKETS ^ ) ; 

if archive then 
begin 
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writeln (outf ilevar) ; 

writeln (outf ilevar /CJ CALCULATION *** + *****^); 
writeln (outf ilevar) ; 
end; 

for i := 1 to act ive.matrix.length do 
begin 

if active.matrix = 1 then 
begin 

max ;= realizat ion.x [i , 1] ; 
for j := 2 to max.intervals do 
begin 

if max < realizat ion_x [i , j] then 
majc := realization.x [i , j] ; 

end; 

cj.matrix [n,max] := cj.matrix [n,m 2 ix] + 1; 
if archive then 
begin 

writeln (outf ilevar) ; 

writeln (outfilevar, 'MAX FOR THIS PASS 
IS ' ,max) ; 

writeln (outf ilevar , 'CUMULATIVE MAX FOR 
I=',max,' IS ' ,cj_matrix[n,max]) ; 

end ; 

end 

else 

begin 

max := realizat ion_y [i , l] ; 
for j:= 2 to max.intervals do 
begin 

if max < realization.y [i , j] then 
max := realizat ion_y [i , j] ; 

end; 

cj .matrix [n, max] := cj.matrix [n,max] + 1; 
if archive then 
begin 

writeln (outfilevar) ; 

writeln (outfilevar, 'MAX FOR THIS PASS 
IS ',maix); 

writeln (outf ilevar , 'CUMULATIVE MAX FOR 
I=',max,' PKTS=',n,' IS 
' , cj .matrix [n ,m2Lx] ) ; 
end; 

end ; 

end ; 
end ; 



{ } 

{ PRINT.CJ } 



{ } 

{ This procedure prints out the matrix of Cj values } 
{ constructed by the procedure GET_Cj . } 

{ } 



^^Ic********* *********************** ♦♦****♦♦♦♦*♦♦♦***♦♦♦♦}- 

procedure PRINT_CJ; 
var 

i, 

j .-integer; 
begin 

writeln I 1234567 89 10 O; 

for i:= 1 to cj.matrix^dim do 
begin 

write (i , ' 10; 

if archive then 

write (outf ilevar , 1 O 10; 

for j:= 1 to cj .matrix. dim do 
begin 

write(^ Ocj_matrix[i , j] O O; 
if archive then 

write (outfilevar, ^ ^ ,cj_matrix [i , j] ) ; 

end; 

writeln; 
if archive then 

writeln (outfilevar) ; 

end ; 

end ; 






{ } 

{ COMB } 

{ } 

{ Given a total number of interferers, every } 

{ combination of interferer packet type must be } 

{ determined. This procedure accepts the total } 

{ number of interferers n and generates every } 

{ combination of the four interferer types which } 

{ add up to the total n. This is basically an } 

{ enumeration of the number of combinations of 4 } 

{ things taken n at a time. This algorithm is } 

{ derived from Reference 7. } 

{ } 



procedure COMB (n ,k ; integer ; var r: index; var mtc ; boolean) ; 
label 10, 11, 15, 20, 30, 31; 
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var 



i, 

j. 
w, 

X, 

y. 

z, 

1 : integer; 
begin 



10 : 



11 : 

15: 



20 : 



30: 



31: 



if mtc = true then 
begin 

goto 20; 

end 

else 

begin 

r[l] := n; 
t := n; 
h := 0; 

for 1 := 2 to k do 
r[l] := 0; 
if r[k] = n then 
mtc := false 
else 

mtc := true; 
goto 31; 
if t > 1 then 
begin 

h := 0; 
end; 

h := h + 1; 
t := r [h] ; 
r[h] := 0; 
r[l] := t - 1; 
r[h+l] := r[h+l] + 1; 
goto 15; 
end ; 

1 := 1 ; 

end; 



{ 

{ GET.COMBINATIONS 

{ 

{ This procedure controls the execution of the 
{ procedure COMB and generates the values w, x, y, 

{ and z. These variables correspond to early-late, 
{ early-early, late-early, and late-late packets 
{ respectively. All combinations given a total 



} 

} 

} 

} 

> 

} 

} 

} 

} 
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{ number of interferers n are placed into the matrix } 

{ vector. } 

{ } 

{:*c:*c;4t:*c:*c*:*c:*c3*c*3*c3*c:*c* + ic:*c3*c3*c:|c5*c:*c54c:*c3tc:*c3*c5*c:tc*************3*c***** + **5*c**}- 

procedure GET_COMBINATIONS ( var vector : comb .matrix; 

n: integer; var vector.length 
integer) ; 



var 

i, 

j> 

w, 

X, 

y> 

z : integer; 
mtc iboolean ; 
v: index; 



begin 

for i := 1 to 4 do 
v[i] := 0; 



{ 



i ;= 0; 
mtc := true; 

if archive then 

writeln (outfilevar, 'COMBINATIONS 

while mtc do 
begin 

if i = 0 then 
mtc := false; 
i := i + 1; 
comb (n,4,v,mtc); 
for j := 1 to 4 do 

vectorEi , j] := v[j] ; 
w := v[l] ; 

X := v[2]; 
y := v[3] ; 
z := v[4] ; 

writeln (w,' ',x,' ',y,' ’ ,z) ;} 
if archive then 
begin 

writeln (outfilevar, w, ' ',x,' ’ ,y , ’ ' ,z 
end; 

end; 

vector_length ;= i; 
if archive then 



writeln 

(outf ilevar , ^******************stc3*c:*cj|cs*c:*c3tc:|cjtcstc;tc:*c*5tcjtc^) ♦ 
end; 



{ } 

{ MAIN } 

{ > 

{ Prograjn execution begins with a query for an } 

{ archive data file. If the archive option is chosen } 

{ all screen I/O will be echoed to the file given. } 

{ The next program query is for the maximum number of } 

{ packets, n. The program will then generate all } 



{ combinations of interferer packet types for a given } 
{ value of n. Each combination is then fed one by } 
{ one into the GENERATE.REALIZATION procedure where } 
{ all possible realizations are generated. This list } 



{ of realizations is then filtered by FILTER to cull } 
{ out invalid realizations. A matrix of valid } 
{ realizations is fed to GET.Cj in order to generate } 
{ the Cj table. This process is repeated for all } 
{ values of n between 1 and the meiximum given by the } 
{ user. } 
{ } 






begin 

archive := false; 

write (^DO YOU WANT A DATA ARCHIVE (y or n) -> 0; 
readln (ans) ; 
if ans = ^y^ then 
begin 

archive := true; 

write (CENTER DATA FILE NAME -> 0 ; 
readln (f ile.name) ; 
assign (outfilevar, f ile_name) ; 
rewrite (outfilevar) ; 
end ; 

write (CENTER THE MAX # OF PKTS PER REALIZATION -> 0 ; 
readln (runs) ; 

initialize.cj .matrix ; 
for n := 1 to runs do 
begin 
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for i := 1 to 220 do 
for j := 1 to 4 do 
vector[i,j] := 0; 

if archive then 
begin 

writeln 

(outf ilevar , ^ 

+ J<e :|c J*c J|C J*C J*C 4: J^C j*e D*cj*e * J*C j*e ^ ) • 

writeln (outf ilevar , ^NUMBER OF PACKETS = \n) ; 
writeln; 
end; 

get.combinations (vector ,n,vector_length) ; 
for i := 1 to vector.length do 
begin 

init.interf erers := vector[i, 1] ^vector [i ,2] ; 
f inal.interf erers := vector [i , 1] +vector[i ,4] ; 
max.interf erers := vector [i,l] + vector [i, 2] + 
vector[i,3] + vector[i,4]; 
min^interf erers := vector[i,l]; 
num. intervals := vector[i,2] + 2*vector [i ,3] + 
vector[i,4] + 1; 

if archive then 
begin 

writeln (outf ilevar) ; 
writeln (outf ilevar) ; 

writeln (outfilevar, ^ 0 ; 

writeln (outfilevar) ; 

writeln (outf ilevar , ^PACKET CONFIGURATION IS \ vector [i , 1] , 

^ \vector[i,2] , ^ ^ , vector [i , 3] , ^ 

^ , vector [i ,4] ) ; 
writeln (outfilevar) ; 

writeln (outf ilevar , ^ INITIAL PKTS = \ init.interf erers) ; 
writeln (outf ilevar , ^FINAL PKTS = ^ ,f inal.interf erers) ; 
writeln (outf ilevar , ^MAX # PKTS = \m 2 LX_interf erers) ; 

writeln (outf ilevar , ^MIN # PKTS = \min.interferers) ; 

writeln (outf ilevar , ^NUM INTERVALS= \num_intervals) ; 
writeln (outfilevar) ; 

writeln (outfilevar, ^ ^ ) ; 

end; 



init_realization_matrix ; 

generate.realizations (init.interf erers , 
max.interf erers , 

num.intervals , min.interf erers , 
active.matrix, 
active_matrix_length) ; 
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filter ( act ive_roatrix. length, 
f inal.interf erers , 

active. matrix) ; 

get.cj (active.matrix, active 
max.interf erers) ; 

end; 

end ; 

print. cj ; 

close (outf ilevar) ; 



num.intervals , 



.rnatr ix.length , 
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APPENDIX B 



DERIVATION OF THROUGHPUT EQUATION 
FOR VARIABLE PACKET ALOHA 

This appendix presents the derivation of the throughput equation for variable 
packet Aloha with power capture. 

As seen in Chapter II equation (2.9), the probability of a tagged packet success- 
fully being transmitted is 

n=l '*• jr=0 

Assign the following identity 

/,(.!) = /(n)2-l"+‘l^i;C,(n) (B.2) 

J=0 

As described in Chapter II equations (2.10) and (2.11), the conditional channel 
departure rate and the channel throughput are given by the following: 

<f){x,u) = g{x + u)ps{x,u) (B.3) 

5 = — / / (f>{x.u)a{x)a{u)dxdu (B.4) 

2 Jo Jo 



Utilizing the binomial identity 



and expanding the exponential term in (B.l), the equation for Ps{x, u) can be rewrit- 



ten as 



1 -t- J2h{n)J2 



n=l k=0 

Expansion of <f){j\u) in terms of pn{x,u) results in 





1 ^k ^ n—k 




1 X U 



(B.6) 



(f>{x, u) = g{x -}- u)e 
Multiplying produces 



1 -f Y, H^) 

n=l k=0 



( ^ ' 


\ k n—k 




H 



(B.7) 



<t>{x.ii) = [gxc-^^c-^'^ + guc-^'^e-^^] 

+ [gxc-^'‘c~^^ + guc~^\-^^]- 

Using this definition of (p{x, u) in the equation for S (B.4) produces two separate 
double integrals. The first is equivalent to the equation (27) in Ref. 2. Assuming 
that Pi(.t) = 0, i> 2 (y) = 0 and that the density function a{x) is exponential, this 
first term for S is found in Ref. 2 to be 



Y Hn) Y 

.n=l k=iO 




(B.8) 



•S'l = 



6’ 



(B.9) 



(^+1)3 

The remaining double integral remains to be evaluated. This term represents the 
increase in throughput due to power capture. The term is written here as 



1 roo roo 

S'2 = - / -h gue-^^e~^^] 

2 Jo Jo 






Ln = l 



k=0 



( " 


1 ^k ^ n — k 




X U 



a{x)a{u)dxdn 



(B.IO) 



Manipulating terms in the above equation produces 
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a[x)a[u)dxdu 



too TOO ^ 

+ / / gue-^^e-^^Yl 

-'0 -/o , „ 



Jc=0 L 



^ 1 
it j 



a(a;)fl(u)(ix(iu > (B.ll) 



Further collection of like terms produces 



S2 = \ k j J • y u"“*’e“®“a(u)du 

+ Z ( /! ) • yT°° 7e-^^a(a:)dx| (B.12) 

Each of the individual integrals above can be solved using the following identity 
[Ref. 5:p. 310]. 



/ .r"c-'^^d.r = n!//-<"+’> (B.13) 

Jo 

Additionally, the following substitutions can be made because packet size is as- 
sumed to be exponentially distributed and because the quantit}' G is given in pack- 
ets/average packet length. 



a(x) = (B.14) 

X 

G = gx (B.15) 

Combining the integral identity in equation (B.13) with the definition of a{x) 
produces this solution for the first integral term in equation (B.12). 




x^+\-9-a{x)dx 



1 r 

X Jo 
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(B.16) 



X X 

_ {k + 1)! + 1 (k+2) 

X X 

(Mjj! (^+2) 

The second term in equation (B.12) is evaluated in the same manner as above. 



/ u^-^e-^'^a{u)du = - / u”-'‘'e"<^+=)Vu 

Jo u Jo 

u u 

^ (»-^’)! ( g^+l )-(n-fc+l) (gi-) 

n u 

In this model, all packets come from the same distribution, therefore x = u. This 
modifies the above equation. 

ir^-^c-^\,{u)du = ^^1^(6' + (B.18) 

The third integral term in equation (B.12) is reduced as follows 

The final integral term is as follows 

r x’^c-^^a{x)dx = (B 20) 

Jo X 

The total throughput of the system is 

.S' = S'l + 52 (B.21) 

Collecting equations (B.9), (B.ll), and (B.16)-(B.20) produces 
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5 = 



(GW S Is, + ‘ > 

•|^(G + l)-'”-‘+'l 



This simplifies to 



(B.22) 



^ 1 oo { n 

S = T^r^ + ^pEMn) E 



(6 + 1)^ - „_J 

(G + 

+ 



"Ht+D * Ti,„ 



(B.23) 



Further simplification produces 



5 = 



G 



1 ■>G C/7 4. 1 '|-(n+3) n 

+ \sT.Hn r^-\L Z ■■'(» + 2 ) 

“ n=l A*=0 



(B.24) 



(G + l)3 z- “ a:- 

The final summation in tlie above equation is over the range of A’ = 0 to oc. 
Because there are no k terms within the summation, the entire summation reduces 
to + 2). Combining terms produces 



5 = 



+ O^' z. TT^TT L -• ' 7 -(n+i) 



(G + 1)^ ■ 



71 1 



(B.25) 



Combining terms and utilizing the identitj" G = gx produces the following solu- 
tion for the throughput of variable packet Aloha with power capture. 



5 = 



G 



N 



(C+l)3 



+ f; + l)-'"+3) ^ c,{n) (B.26) 



71 = 1 



J=0 
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APPENDIX C 



DERIVATION OF THROUGHPUT EQUATION 
FOR SELECTIVE-REPEAT ALOHA 



This appendix presents the derivation of the throughput equation for the variable 
packet Aloha network with pow'er capture and selective-repeat data link control 
presented in Chapter III. 

From Chapter 111 equation (3.2), the probability of tagged packet success to be 



p^{z) = 



3=0 



n=l 



nl 



(C.2T) 



The conditional probability of successful minipacket transmission is given in equa- 
tion (3.3) of Chapter 111 to be 



9('bi/) =P5(-)[1 -Pi(m)][l -P 2 (y)] (C.28) 

Also given in Chapter III is the conditional minipacket departure rate in the interval 
of ?? + 1 minipackets. This equation is rewritten here. 



<?>m(^,y) = ?«y(2 -M)y(^,y) (C.29) 

Equation (3.5) in Chapter III gives the expression for the average minipacket through- 
put and is rewritten here. 



'z 

— _ . r / 4*m{^ 

z + I Jo 



\y)dy 



(C.30) 
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The equation for minipacket throughput given above is now expanded to produce 



the following: 



'z 

Sm = [\ - p 2 {y)]b{y)dy '^{z + l)D{z)e 

z + l Jo ^ 

+=— / [1 - p2{y)]Ky)dy -Y.-Y. fW 

"r z=l n=l 



— mzg 



(n+l) 






n\ 



(C.31) 



]=0 



where D{z) is the probability that a packet consists of z minipackets. 

The first term in the above equation is the same as the equation evaluated in 
Ref. 2 equation (36). The resulting throughput equation for this term is given in 
Ref. 2 equation (37b). Call this part of the total throughput S\. Then from Ref. 2, 



•s. = 



aGv^{2 — r>)[l — pi(m)] 



:{z + l)it'2(l - u)2 



With the following definitions: 



(C.32) 



U' = 1 



V = we ^ 



(C.33) 
(C.34) 

G — mzg (C.35) 

cr = / [I - p 2 {y)]b{y)dy (C.36) 

Jo 

Let S 2 be the remaining term of the total throughput given in equation (C.31). 
By separating terms and noticing that [g[mz + m)]'^ = [gmY{z-]-\Y ^ S 2 is rewritten 



as 



S 2 = 15^11 
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L 

n=l 






N 



/(„)2-C + i>«^ Z ICjin)] J2i- + ir*'D{z)e-""“ 



_ (C.37) 

j=0 z-l 

As stated in Chapter III equation (3.7), assume that D{z) is geometrically dis- 
tributed and is defined as follows: 



^(-^) = ^(1 - (C.38) 

z z 

Utilizing the equations for G, D{z), and w the equation for S 2 can be written as 



follows: 



52 = 



G 



3+1 



[1 — pi(???)]e 



G\n K 






n = l 



77 ! 



(C.39) 



J=0 



2=1 



Assign the following identity: 



l:{„) = (C'.JO) 

^ J=0 

At this point, consolidate the exponential terms and divide out two of the w 
terms to equalize the exponents within the summation over r. 



CO 00 

S 2 — [1 ~ ^ k{n)w~^ X^(- + (C.41) 

^ n=l 2=0 

Assign the following identities; 



k = (. + 1) (C.42) 

t = (n + 1) (C.43) 



This produces the following 
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(C.44) 



^ oo oo 

S 2 = - Pi(”^)]o- Y. Y 

- + ^ n=] k=2 

At this point, modify the summation over k such that the range is from k = 1 
to 00 . To do this, subtract off from the total of the summation the value of the 
summation for k = \ which is v in this case. Therefore, the above equation is 
changed to 



82= ^^[1 -pi{m)]aYHn)w ^ 

^ ^ ^ n=l U=1 



— V 



(C.45) 



Given that |n| < 1, use the following mathematical identity to solve the above 
equation [Ref. 6: p. 142]. 



a. (-■ - <)' 



Using the above identity, the hnal form of 52 is given as 



(C.46) 



52 = 



C 



+ 1 



[1 -pi(m)]CT^ /i(7t) 



-2 



71 = 1 



W 

t r-1 



(1 + 



{r - ky 



— V 



(C.47) 



r=l k=o A'!(/ - 4- + 1)!_ 

The total throughj)ut for selective-repeat Aloha with power capture is given as 



5 = 5i + 52 



(C.48) 
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APPENDIX D 



ANALYTICAL AND NUMERICAL RESULTS FOR 
NEAR/FAR EFFECT CAPTURE 

This appendix provides the details of the analytical and numerical methods used 
to determine the values for Pr{capture\n^ in Chapter III. First, the derivation of 
Pr{capture\n = 1} is provided. This is followed by a MATHCAD file which was 
used to calculate capture for n = 2 and 3. 

A. ANALYTICAL METHOD 

The case of r? = 1 is sufhciently simple to conduct a straightforward mathematical 
derivation of the capture probability given the near/far effect. From Chapter IV the 
following is known: 




(D.49) 



/a'(.t) 



^ -(2/a + l) 

oc;' 



(D.50) 



fviyW) 



^y-(2/- + l) 

lqO 






(D.51) 



For the case of n = 1, both y and x are defined between the values of and 
cxD. Also from Chapter IV, equation (4.10), the definition of fz{z\n). From this 
equation, it is clear that both yz and y must be greater than r““^. given the valid 
ranges for fx{z') and /r(y|?0- This implies that for values of 2 greater than 1, the 
valid range for y is y > Therefore, the equation for fz{^\n) with the proper 



limits of integration is 



fzi = \n^l) 



(D.52) 



/ CO 

yfx{yz)fY{y\n = l)dy 

max 

\aGJ 



= (— 

\aG 



,-(2/a+l) 



ay 



— A/q '' 



Remembering that G = gives us 



(D.53) 

(D.54) 

(D.55) 



gi- 



fz{z\n = 1) 



or; 



_^-(2/a+\) j^-4/a| 



l_-(2/a+l) 



(D.56) 

(D.57) 



a 

To determine the capture probability from the above distribution, simply inte- 
ate from the capture threshold yo to infinity. 



Pr{capture\n = 1} = Pi'{z > 7ob^ = 1} 


(D.58) 


roo 1 


(D.59) 


J'yo ^ 




1 -2/o 

- 2^0 


(D.60) 
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B. NUMERICAL METHODS 



This document is a MATHCAD file used to numerically calculate 
the values of P{capture|n) for values of n from 2 to 3. The 
general approach is to use a numerical convolution to determine 
the n'th convolution of f(y). These convolution results are 
then used in numerical integrations to generate the P{capture|n} 
results. 



Initialize and Define Variables 



M = 1 
(X = 4 



( M corresponds to r ) 

max 



2 

G = T- M 

D=4 (2<D<32 sets the number of convolution points.) 

i s 0 . . 32- D - 1 



2- i + 1 

V = 

i 2D 



f(t) 



T 

2 t 

cx- G 



h(t) = f(t) 




f(t) corresponds to f(w) in the model 



p(t) = f(t) 

The following equation performs a discrete convolution on the 
functions defined as X and Y. 



conv(X,Y) 



1 

D 



■ icfft 



16- D (cfft(X) • cfft(Y) ) 
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step 1. Convolve f(t) with itself to get f(y|n=2) (call this 
C) and plot the resulting functions. 



C s convLh(v) ,p(v) 






Step 2. Convolve f(y|n=2) with f(t) to get f(y|n=3) (called 
D) and plot the relulting function. 





Step 3. Continue in the above manner for as many cases of 
f(y|n) as needed. 
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step 4 . 



Calculate f(z|n) using a linear interpolation of 
the numerical value for f(y|n). 





'2 




'2 


- 


-+1 




30 


- 


ir 


cx 






<x_ 


f (z) := 2 z 


• 


y 




o<- G 




-o< 








A 


2M 





linterp(v,C,y) dy 



Step 5. Calculate the capture probability given n known 
interferers by integrating f(z|n) (f(z) in this 
example) from the capture threshold (3 in this 
example) to a maximum usefull limit which is 
determined by the range of values for which f(y|n) 
is valid. In this case, 30 is the upper limit. 



.30 

Z := f(z) dz 

3 

Z = 0.057 



Step 6. Continue on in the same manner as above to 

calculate the capture probability for different 
values of n. The following equation defines f(z|n) 
for the case where n=3. 





2 




'2 




- 


-+1 




30 


- 




IT 


(X 






oc_ 




g(h) := 2 h 




y 




• linterp(v, D, y) dy 


o<- G 




-cx 








A 


3- M 





. 30 

H := g(h) dh This equation calculates P { capture | n=3 ) 

J 3 

H = 0.022 
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